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Abstract — Recently,  a  3-D  model  (IMPACT35)  was  developed 
to  predict  a  falling  cylindrical  mine’s  location  and  orientation  in 
air-water-sediment  columns.  The  model  contains  the  following 
three  components:  1)  triple  coordinate  transform,  2)  hydrody¬ 
namics  of  falling  rigid  object  in  a  single  medium  (air,  water,  or 
sediment)  and  in  multiple  media  (air-water  and  water-sediment 
interfaces),  and  3)  delta  method  for  sediment  resistance  with  the 
transient  pore  pressure.  Two  mine-impact  burial  experiments 
were  conducted  to  detect  the  mine  trajectory  in  water  column 
[Carderock  Division,  Naval  Surface  Warfare  Center  (NSWC), 
West  Bethesda,  MD,  on  September  10-14,  2001],  and  to  measure 
the  mine  burial  volume  in  sediment  (Baltic  Sea  in  June  2003).  The 
existing  IMPACT35  predicts  a  mine’s  location  and  orientation  in 
the  water  column,  but  not  in  the  sediment  column.  Since  sediment 
resistance  largely  affects  the  mine  burial  depth  and  orientation  in 
sediment,  a  new  method  (bearing  factor)  is  proposed  to  compute 
the  sediment  resistant  force  and  torque.  The  improvement  of 
IMPACT35  with  the  bearing  factor  method  is  verified  using  the 
data  collected  from  the  Baltic  Sea  mine-impact  burial  experiment. 
The  prediction  error  satisfies  near- Gaussian  distribution.  The 
bias  of  the  burial  volume  (in  percent)  prediction  reduces  from 
11%  using  the  delta  method  (old)  to  0.1%  using  the  bearing  factor 
method  (new).  Correspondingly,  the  root-mean-square  error 
(rmse)  reduces  from  26.8%  to  15.8%. 

Index  Terms — Bearing  factor,  burial  depth  and  orientation, 
drag  and  lift  forces  and  torques,  IMPACT35,  mine-impact  burial 
prediction,  sediment  resistance  force  and  torque,  triple  coordinate 
system. 


I.  Introduction 

THE  conclusion  of  the  cold  war  culminated  with  the 
Union  of  Soviet  Socialist  Republics  (U.S.S.R.)  effectively 
ceasing  to  exist  under  international  law  on  December  31,  1991. 
This  historical  event  caused  the  U.S.  military  and  specifically 
the  U.S.  Navy  and  Marine  Corp  Team  to  shift  tactical  emphasis 
from  “blue”  water,  deep-ocean  doctrine  to  littoral  warfare 
doctrine.  This  shift  changed  military  responses  dealing  with  a 
wide  range  of  worldwide  regional  crises  requiring  forward  sea 
basing,  and  expeditionary  force  landing  support. 
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Sea  mines  are  big  threat  in  naval  operations.  Within  the 
past  15  years  three  U.S.  ships,  the  U.S.S.  Samuel  B.  Roberts 
(FFG-58),  Tripoli  (LPH-10),  and  Princeton  (CG-59)  have  fallen 
victim  to  mines.  Total  ship  damage  was  $125  million  while 
the  mines  cost  approximately  $30000  [1].  Mines  have  evolved 
over  the  years  from  the  dumb  “horned”  contact  mines  that 
damaged  the  Tripoli  and  Roberts  to  ones  that  are  relatively  so¬ 
phisticated — nonmagnetic  materials,  irregular  shapes,  anechoic 
coatings,  multiple  sensors,  and  ship  count  routines.  Despite 
their  increased  sophistication,  mines  remain  inexpensive  and 
are  relatively  easy  to  manufacture,  keep,  and  place. 

Accurate  mine  burial  predictions  are  inherently  difficult 
[2],  because  of  unknown  conditions  in  mine  deployment  and 
uncertain  environments  such  as  waves,  currents,  and  sediment 
transports  [3].  The  U.S.  Navy  developed  operational  models 
to  forecast  ocean  environments  for  mine  burial  prediction  [4], 
[5].  Recently,  statistical  methods  such  as  the  Monte  Carlo  [6] 
and  the  expert  system  methods  [3]  have  been  developed.  These 
methods  have  a  core-physical  model  for  falling  rigid  body 
through  air-water-sediment  columns.  The  U.S.  Navy  has  a 
2-D  model  (IMPACT28)  to  predict  a  cylinder’s  trajectory  and 
impact  burial.  The  data  collected  from  the  mine-impact  burial 
experiment  in  the  surf  zone  near  the  Naval  Postgraduate  School, 
Monterey,  CA,  shows  overprediction  of  the  burial  depth  (an 
order  of  magnitude  larger)  using  IMPACT28  [7]. 

A  3-D  model  (IMPACT35)  was  recently  developed  at  the 
Naval  Postgraduate  School  to  predict  a  cylinder’s  trajectory  and 
impact  burial  [8]— [12] .  The  dynamical  system  can  be  simpli¬ 
fied  using  the  following  three  coordinate  systems:  earth-fixed 
coordinate  (E-coordinate),  the  cylinder’s  main-axis-following 
coordinate  (M-coordinate),  and  hydrodynamic  force-following 
coordinate  (F-coordinate).  The  origin  of  both  M-  and  F-coor- 
dinates  is  at  the  cylinder’s  center  of  mass  (COM).  The  body 
forces  and  their  moments  are  easily  calculated  using  the  E-co¬ 
ordinate  system.  The  hydrodynamic  forces  and  their  moments 
are  easily  computed  using  the  F-coordinate.  The  cylinder’s  mo¬ 
ments  of  gyration  are  simply  represented  using  the  M-coordi¬ 
nate.  When  the  mine  penetrates  into  an  interface  between  two 
media  (air-water  or  water-sediment),  the  cylinder  is  decom¬ 
posed  into  two  parts  with  each  one  contacting  one  medium.  The 
body  forces  (such  as  the  buoyancy  force)  and  surface  forces 
(such  as  pressure,  hydrodynamic  force)  are  computed  separately 
for  the  two  parts.  A  fully  3-D  model  is  developed  for  predicting 
the  translation  velocity  and  orientation  of  a  falling  cylindrical 
mine  through  air,  water,  and  sediment.  The  added  value  capa¬ 
bility  of  the  3-D  model  (IMPACT35)  versus  the  2-D  model  (IM- 
PACT28)  is  verified  using  experimental  data. 
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Fig.  1.  M -coordinate  with  the  COM  as  the  origin  X  and  as  the  two 

axes.  Flere.x  is  the  distance  between  theCOV  (B)  andCOM  (X);  ( L.R )  are 
the  cylinder’s  length  and  radius  [8], 


Recently,  two  mine-impact  burial  experiments  were  con¬ 
ducted  to  detect  mine  trajectory  in  the  water  column  [Carde- 
rock  Division,  Naval  Surface  Warfare  Center  (NSWC),  West 
Bethesda,  MD,  on  September  10-14,  2001]  and  to  measure 
the  mine  burial  in  the  sediment  (Baltic  Sea  in  June  2003). 
The  collected  data  are  used  for  model  verification.  Section  II 
describes  basic  physics  of  the  recently  developed  3-D  model 
(IMPACT35).  Section  III  shows  the  added  value  of  I M  PACT35 
in  predicting  mine  movement  in  the  water  column.  However, 
Section  IV  shows  weakness  of  the  existing  IMPACT35  in 
predicting  mine  movement  in  sediment.  Section  V  presents  the 
new  bearing  factor  method  to  compute  the  sediment  resistant 
force  and  torque.  Section  VI  shows  the  improvement  of  the 
bearing  factor  method  in  predicting  mine  burial  in  sediment. 
Section  VII  presents  the  conclusions. 

II.  Description  of  IM  PACT35 

The  3-D  mine-impact  burial  prediction  model  (IMPACT35) 
contains  the  following  major  components:  1)  triple  coordinate 
systems,  2)  momentum  balance,  3)  moment  of  momentum  bal¬ 
ance,  4)  hydrodynamics,  and  5)  sediment  dynamics.  Among 
them,  the  hydrodynamics  (drag  and  lift  forces  and  torques)  have 
been  described  in  [8]  and  [11]— [13],  and  will  not  be  discussed 
here. 

A.  Triple  Coordinate  Systems 

Consider  an  axially  symmetric  cylinder  with  the  center  of 
mass  (COM )  X  and  the  center  of  volume  (COV )  B  on  the  main 
axis  (Fig.  1).  Let  (L,R,x)  represent  the  cylinder’s  length,  ra¬ 
dius,  and  the  distance  between  the  two  points  (X,B).  The  pos¬ 
itive  x-values  refer  to  the  nose-down  case,  i.e.,  the  point  X 
is  lower  than  the  point  B.  Three  coordinate  systems  are  used 
to  model  the  falling  cylinder  through  the  air,  water,  and  sedi¬ 
ment  phases:  earth-fixed  coordinate  (E-coordinate),  main-axis- 
following  coordinate  (M  -coordinate),  and  force-following  coor¬ 
dinate  (F -coordinate)  systems.  All  the  systems  are  3-D,  orthog¬ 
onal,  and  right-handed  [8], 

The  E-coordinate  system  is  represented  by  FE(0,i,j,k) 
with  the  origin  O,  and  three  axes:  x-  and  y- axes  (horizontal) 
with  the  unit  vectors  (i,j)  and  2- axis  (vertical)  with  the  unit 
vector  k  (upward  positive).  The  position  of  the  cylinder  is 
represented  by  the  position  of  the  COM 

X  =  :ri  +  yj  +  zk  (1) 


Fig.  2.  Three  coordinate  systems.  Here,  (i,j,k)  are  the  unit  vectors  of 
E-coordinate  system.  Both  M  -  and  F-coordinate  systems  share  the  same  axis, 
i.e.,  iM  and  iF  are  the  same  unit  vectors  [8]. 


which  describes  translation  of  the  cylinder.  The  translation  ve¬ 
locity  is  given  by 

f  =V,  V  =  (u,v,w).  (2) 

Let  the  orientation  of  the  cylinder’s  main  axis  (pointing 
downward)  be  given  by  iM-  The  angle  between  iM  and  k 
is  denoted  by  tp2  +  n/2.  Projection  of  the  vector  iM  onto 
the  (a?,  2/)-plane  creates  angle  between  the  projection 
and  the  z-axis  (Fig.  2).  The  M -coordinate  system  is  repre¬ 
sented  by  FM(X,iM,jM,kM)  with  the  origin  X,  unit  vectors 
(iM,jM,kM),  and  coordinates  ( xm,vm,zm )■  In  the  plane 
consisting  of  vectors  iM  and  k  (passing  through  the  point  M), 
two  new  unit  vectors  (jM,k M)  are  defined  with  jM  perpendic¬ 
ular  to  the  (iM,k)-plane,  and  kM  perpendicular  to  iM  in  the 
(iM,  k)-plane.  The  unit  vectors  of  the  M  -coordinate  system  are 
given  by  (Fig.  2) 

j M  =  k M  X  ij M,  k M  =  1m  X  jM-  (3) 

The  M  -coordinate  system  is  solely  determined  by  the  orienta¬ 
tion  of  the  cylinder’s  main  axis  iM. 

The  F-coordinate  system  is  represented  by 
FF(X,iF,jF,kF)  with  the  origin  X,  unit  vectors 
(iF,jF,kF),  and  coordinates  (xF,yF,zF).  Let  V„,  be  the 
fluid  velocity.  The  fluid-to-cylinder  velocity  is  represented  by 
Yr  =  VW  -  V,  that  is  decomposed  into  two  parts 

Vr  =  V1+V2,  V1  =  V1iF,  V2  =  V2jF  (4) 

where 

Vi  =  (Vr  •  iF)iF 
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TABLE  I 

Physical  Parameters  of  the  Model  M  ines  in  the  NSWC-Carderock  Experiment  (After  [19]) 


Mine 

Mass 

to) 

p 

(103kg  m'3) 

L 

(m) 

Jl  2 
(kgm2) 

•W3) 

(kg  m2) 

X 

(m) 

1 

16.96 

1.60 

0.505 

0.0647 

0.356 

0 

2 

22.27 

2.10 

0.505 

0.0806 

0.477 

0 

3 

34.93 

1.60 

1.010 

0.1362 

2.900 

0 

4 

45.85 

2.10 

1.010 

0.1696 

3.820 

0 

5 

45.85 

2.10 

1.010 

0.1693 

3.940 

0.0045 

6 

45.85 

2.10 

1.010 

0.1692 

4.570 

-0.077 

Experiment  Impact35 


Fig.  3.  M  ovementof  mine  #6  [L  =  1.01  m,  p  =  2.1  x  103  kg  m~3)  with  \  =  -0.0077  m  and  f2  —  -14.0°  obtained  from  (a)  NSWC-Carderock  experiment, 
(b)  3-D  IM  PACT35  model,  and  (c)  2-D  IM  PACT28  model  (after  [11]). 


is  the  component  paralleling  to  the  cylinder’s  main  axis  (i.e.; 
along  iM),  and 

v2  =  vr  -  (Vr  •  iF)iF 


is  the  component  perpendicular  to  the  cylinder’s  main  axial  di¬ 
rection.  The  unit  vectors  for  the  F-coordinate  are  defined  by 
(column  vectors) 


i/'  =  i  w 


ru 

?"2i  , 

r3i 


jF  =  V2/|V2|,  kF  =  iFxjF.  (5) 


The  F-coordinate  system  is  solely  determined  by  the  orienta¬ 
tion  of  the  cylinder’s  main  axis  (iM)  and  the  water-to-cylinder 
velocity.  Note  that  the  M  -  and  F-coordinate  systems  have  one 
common  unit  vector  iM  (orientation  of  the  cylinder).  Use  of  the 
F-coordinate  system  simplifies  the  calculations  for  the  lift  and 
drag  forces  and  torques  acting  on  the  cylinder. 


[11H13] 


U 

"0" 

V 

=  — 

0 

w 

F  nh  +  F/, 

pn 


(6) 


where  g  is  the  gravitational  acceleration,  n  is  the  cylinder 
volume,  p  is  the  rigid  body  density,  pH  =  m  is  the  cylinder 
mass,  Fnh  is  the  nonhydrodynamic  force  defined  later,  and  Fh 
is  the  hydrodynamic  force  (i.e.,  surface  force  including  drag, 
lift  forces).  Both  Fnh  and  Fh  are  integrated  for  the  cylinder. 
The  drag  and  lift  forces  are  calculated  using  the  drag  and  lift 
laws  with  the  given  water-to-cylinder  velocity  (Vr).  In  the 
F-coordinate,  Vr  is  decomposed  into  along-cylinder  (Vi)  and 
across-cylinder  (V2)  components.  The  nonhydrodynamic  force 
Fnh  is  the  buoyancy  force  (Ft)  for  the  air  and  water  phases 


F  nh  —  Ft  —  k(paII(/,  pwTlg^  (7) 


B.  Momentum  Balance 

The  3-D  translation  velocity  of  the  cylinder  (V)  is  governed 
by  the  momentum  equation  in  the  E-coordinate  system  [8], 


where  ( pa,pw )  are  the  air  and  water  densities  and  Fnh  is  the 
resultant  of  buoyancy  force  (F6),  and  shearing  resistance  force 
(Fs)  for  the  sediment  phase. 
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Experiment  Impact35 


Fig.  4.  M  ovement  of  mine  #5  [L  =  1.01  m,  p  =  2.1  x  103  kg  m-3)  with  \  =  0.0045  m  and  #g  =  42.2°  obtained  from  (a)  NSWC-Carderock  experiment, 
(b)  3-D  IM  PACT35  model,  and  (c)  2-D  IM  PACT28  model. 


Experiment  Impact35 


Fig.  5.  M  ovement  of  mine  #2 [L  —  0.505  m,  p  =  2.1  x  103  kg  m-3)  with  \  =  0  and  =  87.0°  obtained  from  (a)  NSWC-Carderock  experiment,  (b)  3-D 
IMPACT35  model,  and  (c)  2-D  IM  PACT28  model  (after  [11]). 


C.  Moment  of  Momentum  Equation 

The  moment  of  momentum  equation  is  written  in  the  M  -co- 
ordinate  system,  which  rotates  with  the  angular  velocity  of 

Cl  =  1023m  +  u3kM  ■  (8) 

Usually,  the  angular  velocity  around  the  mine’s  main  axis  iM 
(i.e.,  self-spinning  velocity)  is  very  small  (wi  ~  O)  and 


neglected.  Thus,  we  have 

f2  =  cj. 

This  leads  to  zero  centripetal  and  “Coriolis”  terms 

Q  x  (Q  x  u>)  =  0,  — 2J  •(flxu)  =  0.  (9) 
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Fig.  6.  M  odel  verification  from  prediction  of  the  COM  position  x  using  the  NSWC-Carderock  experiment  data  at  several  time  instances:  0.32,  0.64,  0.96,  and 
1.28  s.  H  ere,  the  first  column  is  the  data-l  M  PACT28  comparison,  the  second  column  is  the  data-l  M  PACT35  comparison,  the  third  column  shows  the  histograms 
of  the  model  error  (6x)  for  IM  PACT28,  and  the  fourth  column  shows  the  histograms  of  the  model  error  (6x)  for  IM  PACT35. 


The  inertial  term 


dfi  .  (  dw2  du>3\  .  . 

—  xa,  =  iM(u3-jr-u2-jF)  (10) 

only  has  the  component  along  the  direction  of  iM  (mine’s  main 
axis).  This  term  may  be  neglected  when  the  self-spinning  ve¬ 


locity  is  small.  The  moment  of  momentum  equation  in  the  M  -co- 
ordinate  system  can  be  simplified  by 


J*  —  ~Mnh  +  Mh  (11) 

where  Mn/l  and  are  the  nonhydrodynamic  and  hydrody¬ 
namic  force  torques.  In  the  M  -coordinate  system,  the  moment 
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Scatter  Diagram 
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Fig.  7.  M  odel  verification  from  prediction  of  the  COM  position  .~  using  the  NSWC-Carderock  experiment  data  at  several  time  instances:  0.32,  0.64,  0.96,  and 
1.28  s.  Here,  the  first  column  is  the  data-l  M  PACT28  comparison,  the  second  column  is  thedata-IM  PACT35  comparison,  the  third  column  shows  the  histograms 
of  the  model  error  (6x)  for  IM  PACT28,  and  the  fourth  column  shows  the  histograms  of  the  model  error  (6x)  for  IM  PACT35. 


of  gyration  tensor  for  the  axially  symmetric  cylinder  is  a  diag¬ 
onal  matrix 


The  gravity  force,  passing  the  COM ,  does  not  induce  the  mo¬ 
ment.  The  buoyancy  force  induces  the  moment  in  the  jM  direc¬ 
tion  if  the  COM  does  not  coincide  with  the  COV  (i.e.,  x^O) 


J  = 


■h 

0 

0 


0  0 
J2  o 
0  J3 


(12) 


where  Jlt  J2,  and  J3  are  the  moments  of  inertia.  The  nonhydro¬ 
dynamic  force  usual  ly  contai  ns  the  gravity  and  buoyancy  forces. 


Mb  =  iFftlxcosV’aJM-  (13) 


D.  Sediment  Dynamics 

In  the  existing  IM  PACT 35  model,  the  sediment  resistance  is 
calculated  using  the  delta  method.  This  method  is  based  on  the 
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Fig.  8.  M  odel  verification  from  prediction  of  the  orientation  tp2  using  the  NSWC-Carderock  experiment  data  at  several  time  instances:  0.32,  0.64,  0.96,  and  1.28 
s.  Here,  the  first  column  isthedata-IM  PACT28  comparison,  the  second  column  isthedata-IM  PACT35  comparison,  the  third  column  shows  the  histograms  of  the 
model  error  (<5^2)  for  IM  PACT28,  and  the  fourth  column  shows  the  histograms  of  the  model  error  for  IM  PACT35. 


assumption  that  the  cylinder  pushes  the  sediment  and  leaves 
space  in  the  wake  as  it  impacts  and  penetrates  into  the  sedi¬ 
ment.  This  space  is  refilled  by  water  and  the  water  cavity  is  pro¬ 
duced.  At  the  instant  of  penetration,  the  total  resistant  force  on 
the  cylinder  is  represented  by  [16]— [18] 


namic  forces  (per  unit  area)  at  the  point  r  over  the  cylinder’s 
surface,  crsed  is  the  area  of  the  cylinder’s  surface  below  the 
water-sediment  interface,  Fpw  is  the  pore  water  pressure  force 
on  the  whole  cylinder,  and  ^-function  is  defined  by 


F s  =  J  [sm  +  fsh)  +  f?  +  f%]da  +  Fpw  (14) 

^Tsed 

where  (f6s,f sh)  and  are  the  sediment  buoyancy  and 

shear  resistance  forces  and  water  buoyancy  and  hydrody¬ 


r  1  v  •  n  >  0 
10  v  •  n  <  0 


(15) 


which  shows  that  the  sediment  buoyancy  and  shear  resistance 
forces  act  when  the  cylinder  moves  towards  them.  Here,  v  is 
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Fig.  9.  Temporal  evolution  of  model  performance  evaluated  using  the  NSWC-Carderock  data:  (a)  observational  data  number,  (b)  rmseof  s,  (c)  rmseof  zt  (d)  rmse 
of  ip2,  (e)  STD  of  6x,  (f)  STD  of  6z,  and  (g)  STD  of  with  the  solid  curves  for  I M  PACT28  and  dashed  curves  for  I M  PACT35  from  (b)  to  (g). 


the  velocity  at  point  r  (represented  in  the  M  -coordinate)  on  the 
cylinder  surface 

v  =  V  +  ui  x  r.  (16) 

III.  Verification  of  IM  PACT35  in  the  Water  Column 

The  NSWC-Carderock  experiment  was  conducted  on 
September  10-14,  2001  in  the  Explosion  Test  Pond,  which  is 


the  only  explosive-related  test  pond  in  the  United  States  with 
the  capability  of  providing  high-speed  underwater  photography 
given  its  exceptional  water  clarity.  In  addition,  the  facility’s 
concrete  floor  thickness  and  reinforcement  is  sufficient  to  allow 
impact  of  45-kg  cylinders  without  additional  floor  protection. 
The  pond  in  plan  view  is  a  regular  pentagon  with  each  side  of 
41  m.  During  the  experiment,  six  model  mines  (Table  I)  with 
mass  varying  from  16.96  to  45.85  kg  were  released  to  the  pond 
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Fig.  10.  Optical  mine  used  in  the  Baltic  Sea  experiment  (after  [24]). 


TABLE  II 

Physical  Parameters  of  the  Full-Size  Optical  Mine  in  the 
Baltic  Sea  Experiment  (After  [24]) 


Length  (m) 

1.466 

Diameter  (m) 

0.470 

Taper  Diameter  (m) 

0.395 

Taper  Length  (m) 

0.150 

Volume  (m3) 

0.252 

Mass  (kg) 

550.00 

Z  (m) 

0 

j,(W) 

14.82 

h,  fitam2) 

105.00 

TABLE  III 

' sica l  Parameters  of  the  E 

iNVIRONMENT  IN 

Baltic  Sea  Experiment 

(After  [24]) 

Air  Density  (kg  m'3) 

1.22 

Water  Density  (kg  m'3) 

Air  Kinetic  Viscosity  (m2  s"1) 

1025.8 
1.46x1  O'5 

Water  Kinetic  Viscosity  (m2  s 

T  1 . 1 3  x  1 0"6 

with  the  water  depth  at  7.92  m  [19],  [20].  The  data  set  collected 
from  the  NSWC-Carderock  experiment  was  used  to  evaluate 
the  added  value  of  I M  PACT35  versus  I M  PACT28  (2-D  model). 

A.  Near  Horizontal  Release 

M  odel  mine  #6  was  released  to  the  water  with  02  =  -14° 
(near  horizontal,  see  Fig.  2).  The  physical  parameters  of  this 
mine  are  given  by 

L=1.01m,  p=2.10xl03  kg  m-3,  x=-°'077m, 

rn=45.85  kg,  Ji  =  0.1692  kg  m2,  J2  =  J3  =  4.570  kg  m2. 

(17) 

The  initial  conditions  are  given  by 

x0=yo  =  zo  =  O,  u0  =  v0  =  w0  =  0, 
r’V’io  =  0,  "020  =  —14°,  t/>3o  =  0,  u>io  =  w2o  =  W30  =  0. 

(18) 


Substitution  of  the  model  parameters  (17)  and  the  initial 
conditions  (18)  into  IM  PACT28  and  IM  PACT35  leads  to  the 
prediction  of  the  mine’s  translation  and  orientation  that  are 
compared  with  the  data  collected  during  the  experiment  at  each 
time  step  (Fig.  3).  The  new  3-D  model  (IM  PACT35)  simulated 
trajectory  agrees  well  with  the  observed  trajectory.  Both  show 
the  same  pattern  and  the  same  travel  time  (1.92  s)  for  the 
cylinder  passing  through  the  water  column.  Flowever,  the  3-D 
model  (I  MPACT35)  is  better  than  the  2-D  model  (IMPACT28) 
in  predicting  the  mine’s  movement  in  the  water  column. 

B.  Near  45°  Release 

M  odel  mine  #6  was  released  to  the  water  with  02  =  42.2°. 
The  initial  conditions  are  given  by 

xo=yo  =  zo  =  ®,  uo=v0  =  w0=0, 

V;io=0,  '020  =  42.2°,  030  =  0,  wio  =  w2o  =  W3o  =  0.(19) 

Substitution  of  the  model  parameters  (17)  and  the  initial 
conditions  (19)  into  IMPACT28  and  IMPACT35  leads  to  the 
prediction  of  the  mine’s  translation  and  orientation  that  are 
compared  with  the  data  collected  during  the  experiment  at  time 
steps  (Fig.  4).  Both  3-D  model  (IMPACT35)  and  2-D  model 
(I  M  PACT28)  simulated  trajectories  and  travel  times  agree  well 
with  the  observed  trajectory. 

C.  Near  Vertical  Release 

Model  mine  #2  was  released  to  the  water  with  02  =  87°.  The 
physical  parameters  of  this  mine  are  given  by 

L  =  0.505  m,  p  =  2.10  x  103  kg  m-3,  x  =  0, 
to  =  22.27  kg,  Ji  =  0.0806  kg  m2, 

J2  =  J3  =  0.477  kg  m2.  (20) 

The  initial  conditions  are  given  by 


xo=yo  =  z0  =  0,  u0  =  vo  =  wo=0, 

'i/’io=0,  020  =  87°,  030  =  0,  wio  =  o;2o  =  W3o  =  0.  (21) 

The  predicted  cylinder’s  translation  and  orientation  are  com¬ 
pared  with  the  data  collected  at  time  steps  (Fig.  5).  The  3-D 
model  (IM  PACT35)  simulated  trajectory  agrees  well  with  the 
observed  trajectory.  Both  show  the  same  straight  pattern  and 
the  same  travel  time  (1.84  s)  for  the  cylinder  passing  through 
the  water  column. 

D.  Statistical  Error  Analysis 

Figs.  6-8  show  scatter  diagrams  and  histograms  for  pre¬ 
dicting  the  mine’s  location  and  orientation  [x(t), 

In  the  scatter  diagrams,  the  points  cluster  around  the  diagonal 
line  of  <zm  =  <;0  using  IMPACT35  (second  column),  and  the 
points  are  spreading  out  of  the  diagonal  line  using  IM  PACT28 
(first  column),  which  confirms  that  IM  PACT35  predicts  COM 
position  more  accurately  than  IMPACT28.  Histograms  of 
model  errorsforCOM  position  have Gaussian-type distribution 
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Density  (kg  m  3) 

Fig.  11.  Sediment  density  pj  )  and  static  shear  strength  5(s)  profiles  in  the  Baltic  Sea  from  the  cores  collected  at  59  mine-impact  sites  during  the  mine-drop 
experiment  in  J  une  2003:  (a)  individual  density  profiles,  (b)  individual  static  shear  strength  profiles,  (c)  mean  density  profile,  and  (d)  mean  static  shear  strength 
profile  (after  [24]). 


with  near-zero  mean  and  small  standard  deviation  (STD)  using 
IMPACT35  (fourth  column),  and  non-Gaussian-type  distribu¬ 
tions  with  large  STD  using  IM  PACT28  (third  column). 

The  total  number  of  observational  points  at  each  time  sample 
N(t)  is  around  41  as  t  <1.2  s  and  reduces  quickly  with  time 
as  t  >1.2  s  Fig.  9(a),  which  indicates  that  the  model  verifica¬ 
tion  is  reliable  for  t  <1.2  s.  The  added  value  of  IM  PACT35 
is  easily  seen  from  Fig.  9.  For  example,  the  root-mean-square 
error  (rmse)  of  COM  prediction  is  much  smaller  when  using 
IM  PACT35  than  IM  PACT28  [Fig.  9(b)  and  (c)].  The  STDs  of 
the  model  errors  for  the  COM  prediction  are  also  much  lower 
when  using  IMPACT35  than  IMPACT28  [Fig.  9(e)  and  (f)]. 
The  rmse  of  orientation  prediction  is  smaller  when  using  IM- 
PACT35  than  IM  PACT28  for  t  <1  s.  When  t  >1  s,  the  rmse 
of  ip 2  is  around  half  when  using  IM  PACT35  than  when  using 
IMPACT28  [Fig.  9(d)  and  (g)]. 


IV.  Verification  of  IM  PACT35  in  Sediment 

The  Baltic  Sea  experiment  was  conducted  in  June  2003  by 
the  German  Federal  Armed  Forces  Underwater  Acoustic  and 
M  arine  Geophysics  Research  Institute  (FWG,  Kiel,  Germany) 
[21]  with  the  full -size  optical  mine  (Fig.  10)  which  is  allowed  to 
free  fall  from  the  winch.  Table  1 1  shows  the  physical  parameters 
of  the  optical  mine.  Table  III  lists  the  physical  environments  in 
the  Baltic  Sea.  The  full-size  optical  mine  was  released  59  times. 
T  he  water  depths  of  the  drop  sites  were  between  25.0  and  26.5  m 
(Fig.  11).  The  volume  of  mine  burial  percentage  was  measured. 
The  readers  are  referred  to  [6]  for  detailed  information. 

After  running  IM  PACT35  with  the  delta  method  for  the  sedi¬ 
ment  resistance  (14)  for  each  gravity  core  regime  [ps(z),  S(z)], 
the  predicted  and  observed  burial  volumes  (in  percent)  were 
compared  (Fig.  12).  The  bias  (mean  predicted  minus  observed 
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Fig.  12.  Scatter  diagrams  of  buried  volume  prediction  against  the  Baltic 
Sea  experimental  data.  Here,  IM  PACT 35  uses  the  delta  method  for  sediment 
resistance. 


values)  and  rmse  of  burial  volume  are  11%  and  26.8%.  The 
correlation  coefficient  between  predicted  and  observed  burial 
volumes  is  0.374%. 

The  histograms  of  the  burial  volume  (in  percent)  are 
very  different  between  the  Baltic  Sea  experiment  and  the 
model  prediction  of  IMPACT35  using  the  delta  method 
[Fig.  13(a)  and  (b)].  In  the  experiment,  the  probability  density 
function  (pdf)  has  a  peak  at  burial  volume  of  50%  with  a 
frequency  of  26.  However,  the  model  predicted  pdf  has  a  peak 
at  30%  with  a  frequency  of  11. 

V.  Bearing  Factor  M  ethod 

A  new  (bearing  factor)  method  is  presented  for  calculating 
the  sediment  resistant  force  and  torque. 


Buried  Volume  (%) 

(b) 

Fig.  13.  Histograms  of  buried  volume  (in  percent)  from  (a)  Baltic  Sea 
experiment  and  (b)  prediction  using  IM  PACT35  with  the  delta  method. 


(outward  positive),  and  zws  represents  the  vertical  coordinate  of 
the  water-sediment  interface. 

The  shear  resistance  force  (Ff )  is  in  the  opposite  direction 
of  v  and  acts  on  the  mine.  Its  magnitude  is  proportional  to  the 
product  of  the  sediment  shear  strength  ( S )  and  the  rupture  area 
(A,  projection  of  sediment-contacting  area  perpendicular  to  the 
velocity  V)  with  a  nonnegative  bearing  factor  N  [15] 


A.  Sediment  Resistance 


F*  =  —tN(p,v)SA,  N>  0.  (24) 


When  the  mine  impacts  and  penetrates  into  the  sediment,  it 
creates  a  I  arge  transi  ent  pore  pressure  i  n  the  sedi  ment  that  causes 
ruptures  in  the  sediment  and  influences  the  resistance  force  on 
the  cylinder  [22],  [23],  The  resistance  of  the  sediment  to  the 
mine’s  penetration  isassumedto  have  the  foil  owing  three  com¬ 
ponents:  buoyancy  force  Ff,  hydrodynamic  (drag  and  lift)  force 
Fsh  (similar  to  air  and  water),  and  shear  resistance  force  Ff  (re¬ 
sistance  to  the  rupture) 

Fs  =  Fsb  +Fsh+Fsr.  (22) 


The  sediment  resistance  torque  includes  the  hydrodynamic  and 
shearing  resistance  torques 

Ms  =  r  x  (Fj!  +  Fsh)  =  -  [N(p,v)SA]r  xr  +  r  xF|.  (25) 

Here,  p  is  the  nondimensional  penetration  depth  scaled  by  the 
diameter  of  the  cylinder  (2R).  The  sediment  density  in  (23)  and 
shear  strength  ( S )  in  (24)  and  (25)  are  measured. 

The  bearing  factor  increases  with  p  and  decreases  with  the 
decreasing  speed 


The  sediment  buoyancy  force  per  unit  area  is  defined  by 

zws 

F|  =  -n  /  ps(z')gdz'  (23) 

Z 

where  ps(z)  isthesedimentwet density  (usually  obtained  from 
the  sedi  ment  data),  n  is  a  unit  vector  normal  to  the  mine  surface 


N(p,v)  =  [pi  y/12] 


1  +  A  log 


(26) 


where  A  is  the  ^-effect  parameter,  (g\  .//,2)  are  the  p- effect  pa¬ 
rameters  [15],  and  vcii  is  the  critical  speed.  The  bearing  factor 
N  is  amplified  [1  +  Alog^/^cn)  >  1]  for  v  >  vcri,  indepen¬ 
dent  of  v  for  v  =  vcri,  and  is  reduced  [1  +  Alog(t>/t)cri)  <  1] 
for  v  <  vcii.  During  the  penetration,  v  decreases  since  the 
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Fig.  14.  Calculation  of  immersed  area  (after  [12]). 


buoyancy,  hydrodynamic,  and  shear  resistance  forces  oppose 
the  penetration.  Decrease  of  v  reduces  bearing  factor  N.  The 
two  p-effect  parameters  are  given  by  [15] 

Hi  =  9.3  H2  =  0.7. 

N  ote  that  the  nondimensional  (26)  and  the  two  p-effect  param¬ 
eters  are  derived  from  a  small  probe  under  axial  impact  condi¬ 


tion.Validity  of  the  bearing  factor  method  to  mine-impact  burial 
needs  thorough  evaluation. 

Since  N  cannot  be  negative,  when  v  decreases  to  a  critical 

V  =  w~(1/A)  (27) 

the  bearing  factor  N(p,v )  and  in  turn  the  shearing  resistance 
force  become  zero.  T  he  mi  ne  ceases  the  penetrati  on  i  n  sedi  ment. 
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N  ote  that  A  and  vcri  are  the  two  tuni  ng  parameters  of  the  numer- 
ical  model.  In  this  paper,  we  use 

A  =  0.2  vCTi  =  0.0015  m  s_1.  (28) 

B.  Rupture  Area 

Usually,  after  passing  through  the  water  column,  the 
cylinder’s  velocity  reduces.  The  rupture  area  can  be  repre¬ 
sented  approximately  by  the  contact  surface  area.  Let  ( Af ,  At) 
be  the  sediment-contacting  areas  on  the  circular  bottom  (fan 
shape)  and  side  (curved  trapezoid).  The  rupture  area  A  is 
the  summation  of  the  projections  of  (Af,At)  on  the  plane 
perpendicular  to  the  mine’s  velocity  V 


A  =  A\  +  A2  A\  =  Af  cos  a  ^2  =  sin  a:  (29) 

where  a  is  the  angle  between  V  and  the  mine’s  axis.  Let  R  be 
the  radius  of  the  cylinder,  h  be  the  depth  of  Af,  and  (hi,h2) 
be  the  depths  of  the  trapezoid  At  immersed  in  the  sediment 
(Fig.  14).  Letthelength  of  the  trapezoid  b  eLt.  Af  is  computed 
simply  by 


Fig.  15.  Scatter  diagrams  of  buried  volume  prediction  against  the  Baltic  Sea 
Experiment  data.  Here,  IM  PACT35  uses  the  bearing  factor  method  for  sediment 
resistance. 


Af  =  6-R2-(R-h) 0  =  cos-1  f  ■ 

(30) 

At  is  computed  by 


2LtR. 


for  hi  >  h2  >  R 


- rrpm-02 ) 


(hi-h2y 

At={  -(R-h\)bi  +  (R-h2)b2],  for  h2<hi<R 

2LtlR+-^— 

( R~h2 ) 

„  x  [R2  (f-02)  +  (R-h2)b2l  ,  for  hi  >R>h2 

(31) 

where  (Ltl,Lt2)  are  the  lengths  of  subcylinders  with  (h  > 
R,  h  <  R)  and  61  and  b2  are  computed  by 


h  =  b2  =  y/R?  -  (R  -  hi)*.  (32) 


Buried  Volume  (%) 

Fig.  16.  Histogram  of  buried  volume  (in  percent)  prediction  using  IM  PACT35 
with  the  delta  method. 


VI.  Model  Improvement  Using  the  Bearing 
Factor  M  ethod 

M  odel  improvement  is  evaluated  using  the  Baltic  Sea  exper¬ 
iment  data.  Similar  to  the  process  described  in  Section  V,  the 
modeled  (with  the  bearing  factor  method)  and  observed  burial 
volumes  (in  percent)  were  compared  (Fig.  15).  As  evident,  the 
sediment  resistance  using  the  bearing  factor  method  improves 
the  prediction  capability  of  IMPACT35.  The  bias  (mean  pre¬ 
dicted  minus  observed  values)  and  rmse  of  the  burial  volume 
reduce  to  0.1%  and  15.8%.  The  correlation  coefficient  between 
observed  and  predicted  burial  volumes  increases  to  0.435  (using 
the  bearing  factor  method)  from  0.374  (using  the  delta  method). 

The  histogram  of  the  predicted  burial  volume  (in  percent) 
using  the  bearing  factor  method  (Fig.  16)  is  closer  to  that  of 


the  Baltic  Sea  experiment  [Fig.  13(a)]  than  the  histogram  of 
the  predicted  burial  volume  (in  percent)  using  the  delta  method 
Fig.  13(b).  In  the  experiment,  the  histogram  hasapeak  at  burial 
volume  of  50%  with  a  frequency  of  26.  In  the  prediction,  the 
histogram  has  a  peak  at  30%  with  a  frequency  of  11  using  the 
delta  method  Fig.  13(b).  However,  the  histogram  has  a  peak  at 
50%  with  a  frequency  of  17  (Fig.  16). 

The  histogram  of  the  model  error  is  more  symmetric  around 
the  zero  error  using  the  bearing  factor  method  than  using  the 
delta  method.  For  example,  the  histogram  has  a  peak  at  zero 
error  with  frequency  of  16  using  the  bearing  factor  method 
[Fig.  17(b)],  and  a  peak  at  -10%  error  with  the  frequency 
of  12  using  the  delta  method  [Fig.  17(a)].  The  burial  volume 
prediction  is  unbiased  using  the  bearing  factor  method,  and 
negatively  biased  using  the  delta  method. 
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(a) 


Buried  Volume  (%) 

(b) 

Fig.  17.  Histograms  of  the  buried  volume  prediction  error  of  IM  PACT35  using 
(a)  the  delta  method,  and  (b)  the  bearing  factor  method. 


VII.  Summary 

The  following  points  summarize  this  paper. 

1)  The  3-D  mine-impact  burial  prediction  model 
(IMPACT35)  was  recently  developed  to  predict  the 
translation  and  orientation  of  falling  cylindrical  mine 
through  air,  water,  and  sediment.  It  contains  the  foil  owing 
three  components:  triple  coordinate  transform,  cylinder 
decomposition,  and  hydrodynamics  of  falling  rigid  object 
in  a  single  medium  (air,  water,  or  sediment)  and  in 
m  u  I  ti  p  I  e  m  ed  i  a  ( ai  r-w  ater  and  w  ater-sed  i  men  t  i  n  terf  aces) . 

2)  Data  collected  from  two  mine-impact  burial  experiments 
were  used  to  verify  the  existing  IMPACT35.  The  pre¬ 
dicted  mine  track  and  orientation  in  the  water  column 
agree  quite  well  with  the  N  SW  C-Caderock  data.  The  rmse 
of  the  mine’s  position  and  orientation  is  much  smaller 
using  IMPACT35  than  using  IMPACT28.  However,  the 
predicted  mine  burial  volume  in  the  sediment  was  not  as 
good  as  the  mine  trajectory  in  the  water  column. 

3)  Calculation  of  the  sediment  resistant  force  and  torque  is 
updated  from  the  delta  method  (old)  to  the  bearing  factor 
method  (new).  With  the  bearing  factor  method,  the  predic¬ 
tion  capability  of  IM  PACT35  has  been  greatly  improved. 
The  prediction  error  satisfies  near-Gaussian  distribution. 
The  bias  of  the  burial  volume  (in  percent)  prediction  re¬ 
duces  from  11%  using  the  delta  method  (old)  to  0.1%  using 


the  bearing  factor  method.  Correspondingly,  the  rmse  re¬ 
duces  from  26.8%  to  15.8%. 

4)  I M  PACT35  developed  in  this  paper  is  for  cylindrical  mines 
only.  It  is  necessary  to  extend  the  modeling  effort  to  more 
realistic  mine  shapes  such  as  Rockan  and  M  anta  for  oper¬ 
ational  use. 
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